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ABSTRACT 



Context. There is general agreement that supemovae la correspond to the thermonuclear runaway of a white dwarf that is part of 
a compact binary, but the details of the progenitor systems are still unknown and much debated. One of the proposed progenitor 
theories is the single-degenerate channel in which a white dwarf accretes from a companion, grows in mass, reaches a critical mass 
limit, and is then consumed after thermonuclear runaway sets in. However, there are major disagreements about the theoretical delay 
time distribution and the corresponding time-integrated supernova la rate from this channel. 

Aims. We investigate whether the differences are due to the uncertainty in the common envelope phase and the fraction of transferred 
mass that is retained by the white dwarf. This so-called retention efficiency may have a strong influence on the final amount and timing 
of supemovae la. 

Methods. Using the population synthesis code SeBa, we simulated large numbers of binaries for various assumptions on common 
envelopes and retention efficiencies. We compare the resulting supernova la rates and delay time distributions with each other and 
with those from the literature, including observational data. 

Results. For the three assumed retention efficiencies, the integrated rate varies by a factor 3-4 to even more than a factor 100, so in 
extreme cases, the retention efficiency strongly suppresses the single-degenerate channel. Our different assumptions for the common 
envelope phase change the integrated rate by a factor 2-3. Although our results do recover the trend in the theoretical predictions from 
different binary population synthesis codes, they do not fully explain the large disagreement among them. 

Key words, binaries: close, symbiotic ~ white dwarfs - supemovae: general - novae 



1. Introduction 

Supernova la (SNIa) light curves are scalable to one proto- 
type light curve due to th e consistent pr oduction of a cer- 
tain peak luminosity (e.g. iPhillipsI Il993h . This characteris- 
tic makes it relatively easy to estimate distances to observed 
SNela, which is why th ey can be used as st a ndard can- 
dles i n cosmology (e . g. Leibundgut et al.l 119911: iRiess et al] 
Il998t iPerlmutteretall Il999l) . SNela also stronglv affect the 



Galactic chemical evolution throu gh the expulsion of iron (e.g. 
Ivan den Bergh & Tammanrilll991h . In addition, SNela play an 
important role in astroparticle physics, because the accompany- 
ing shoc ks are prime accelerator site s for galactic cosmic ray 
■ particles dBlandford & Ostrikeill978l ). Although they are part of 
many research fields, their origin and the details of the underly- 
ing physical processes are not fully understood. 

It is generally accepted that SNIa events are caused by the 
thermonuclear explosions of carbon/oxygen (C/O) whit e dwarfs 
(WDs) with masses near the Chandrasekhar mass (iNomotol 
Il982ah . Of the two classical progenit or scenarios, we look 
at the single-degenerate (SD) channel (IWhelan & IbenI [T973h . 
in which a WD accretes from a companion. The double- 
dege nerate (DD) channel describes the merger of two WD s 
(e.g. llben & Tutukovlll984l:IWebbinklll984l:lToonen e"tani2012l) . 
Other channels are amongst others accreting WDs from hehum- 



rich, non-degenerate companions (see e.g. I Wang et ani2009bl) . 
and merger s between a WD and t he core of an asymptotic giant 
branch star (lKashi&Sokeiil2011l) . 

The progenitor theory should give a reliable description of 
the evolution of a SNIa progenitor binary. Beyond this, it should 
be able to explain and reproduce general features of the type la 
supernova class. Both the SD and DD scenarios have problems 
in matching observed features of the SNela even ts. We briefly 
men tion the most impor tant issues; however, see iLiviol (l2000h 
and IWang & Hani (I20T2I) for reviews. Regarding the DD sce- 
nario, a serious concern is whether the collapse of the remnant 
would lead to a superno va or to a neutron star through accretion- 
induced collapse (see 'Nomoto & IbenI 19851; ISaio & Nomotol 
1985: Versanti et al. 2003; Yoon et al.l2()ol Pakmor et a l.l2010l 
2012ilShen et al.ll2012l) . In the SD channel, a SNIa-like event 
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is more easily reproduced in the simulations of the explosion 
process, although the explosion process needs to be fine-tuned 
to reproduce the observed spectra and lightcurves (see e.g. 
^Hillebrandt & Niemever 2000i, for a review). Another issue is 
the long phase of supersoft X-ray emission SNIa SD progeni- 
tors should go through. It is unclear whether there are enough 
of these sources to account for the SNIa rate (see Di Stefanol 
120101: iGilfanov & Bogdanll2010l;lHachisu et al.ll2010l) . Archival 
data of known SNela have not shown this e mission unambigu- 
ously , although there is pos s ibly one case ( see \^ss & NelemansI 
I2008t iRoelofs et al.1 120081: iNielsen etalJ l2010l) . Finally, sev- 
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eral predictions for the SD rate have been made by differ- 
ent groups using binary population synthesis (BPS) simula- 
tions (fVung elson & Liv io 200 0; Han & Podsiadlowski 2004 



Ruiter et aO 12009: M eng et all | 2009: Mennekens et al. I201L 
Wang et al.ll2010^ .YungelsonI 12010; .Claevs et al.. ,201 L) . time- 



integrated rates, as well as the distribution of rates over time, 
the so-called del ay t ime d istribution (DTD). For an overview, 
see iNelemans et al] (|2012|) . However, the results show a wide 
spread and do not agree with each other or with observational 
data. The exact origin of these differences is so far unclear. 

In this paper we try to uncover the reason for these differ- 
ences, and we focus in particular on the efficiency with which 
the WD accretes and retains matter from its companion and on 
the co mmon envelope (CE) phase. The re are several prescrip- 
tions jNomoto et al.ll2007l: iRuiteret al.1 [20091: lYungelsonl 12010. 
hereafter NSKH07, RBF09 and YIO respectively) in the litera- 
ture for the efficienc y of retairiing matt er by the WD. The re- 
tention efficiencies of lNSKHOTllRBFOgl and lYld differ strongly. 
We study how the retention efficiency influences the SNIa rate 
and the DTD. Furthermore we investigate if the differences in 
the assumed retention efficiencies can explain the differences in 
the predicted SNIa rate from different BPS codes. 

Differences between the retention efficiencies arise from the 
uncertainty in the novae phase and the strengt h of the optically 
thick WD wind at high mass transfer rates (iKato & Hachisul 
Il994l) . The optically thick wind s tabilises mass transfe r such 
that a CE phase can be avoided. (iHachisu et al.l [T999bl here- 
after HKN99b) argues that the WD wind can interact with the 
envelope of the companion and strip some of the envelope 
mass from the donor, which stabilises the mass transfer further. 
The wind-stripping effect affects the retention efficiency mod- 
erately, but enlarges the parameter space for producing SNela 
effec tively. However, for low metallicities the wind atten uates 
(e.g. iKobavashi et"al]|1998al: iKobavashi & Nomotdl2009h . but 
a low-metaUicity threshold of SNela in comparison with SNe 
type II have not bee n found in observations (Prieto et al. 200^ 
iBadenes et ai]|2009l) . Regarding the novae phase, despite much 
progress in the understanding of classical novae, they are still 
poorly unders tood e.g. the mixing betw een the accreted envelope 
and WD (e.g. lDenissenkov et al.ll2013h . Currently it is unclear if 
a cycle of novae outbursts removes more mass from the WD 
than the accreted envelope, effectively reducing the WD in mass 
dPrialriikl [19861 iPrialnik & Kovetzl [l99l iTownslev & BildsterJ 
12004 lYaron et al. I l2005h . In particular, the results from helium 
novae are different, which affects the chain of nuclear burning in 
the SD channel for which first hydrogen-rich material is burned 
into helium-rich material and consequently the helium-rich mat- 
ter into carbon-rich material. The uncertainty in helium accretion 
creates the strongest difference between the assumed retention 
efficiencies. 

In Sect.[2] the general evolution of a progenitor binary wifl 
be described. The CE phase and the growth in mass of the WD 
will be discussed in detail. To predict the distribution of SNIa 
rates for different assumptions we simulate the evolution for a 
large population of binaries using the population synthesis code 
SeBa, which is described in Sect.[3] Different methods for select- 
ing those binaries that will evolve into a SNIa are also outlined. 
This includes a straightforward selection on binary parameters 
and an implementation in SeBa of the efficiency of mass reten- 
tion on the WD. The resulting rates are presented in Sect.[4] and 
compared with theoretical predictions from different groups in 
Sect.[5] 



2. Evolution of single degenerate SNIa progenitors 

2.1. Global evolution 

The evolution of most binaries starts with two zero-age main se- 
quence (ZAMS) stars in orbit around their common centre of 
mass. This is shown as the first stage in Fig.[T] which shows the 
evolutionary stages a SD progenitor binary goes through. The 
initially more massive (hereafter primary) star has a mass of 
about 3-8M0 and evolves into a C/O WD of about 0.6-1.2Mo. 
In order to form a compact system, binaries go through one or 
several mass transfer phases. The mass transfer can be stable or 
unstable, where the latter is described by the CE phase, the onset 
of which is shown in the second stage in Fig.[T] Here the primary 
has evolved off the main sequence (MS) and overfills its Roche 
lobe. In the following CE phase the envelope of the primary star 
engulfs the initially less massive star, hereafter secondary. The 
core of the primary and secondary spiral inward through the en- 
velope and the envelope itself is expelled. After the CE phase 
the binary consists of a degenerate C/O WD and a MS star in a 
close binary, see the third stage in Fig.[T] Note that the evolution- 
ary stages following the MS are relatively short, so that even an 
initially small mass difference can lead to a binary with a WD 
and a MS star. In spite of the importance of the CE phase for 
creating short period systems containing compact objects, the 
phenomenon is not yet well understood. In Sect. l2.2l we discuss 
two prescriptions for the CE phase. 

If the orbital separation of the binary after the CE phase is 
small, the secondary will fill its Roche lobe when it evolves. 
Depending on the exact orbital distance this may occur when 
the secondary is a late MS star, a Hertzsprung gap star, or after 
it has evolved into a red giant. 

Mass from the outer layers of the secondary will be trans- 
ferred through the first Lagrangian point onto the WD. The 
hydrogen-rich material accumulates in a layer on the surface of 
the WD. When the pressure and temperature are high enough, 
the hydrogen layer ignites. As a result (a fraction of the) mass 
can be expelled from the WD and (a fraction of the) mass can 
be retained by the WD to form a helium layer on the WD. If the 
helium layer ignites, carbon and oxygen will be formed. Under 
the right conditions this mass too can be retained by the WD. 
A favourable combination of these two stages is needed for the 
C/O WD to be able to grow significantly in mass. In section [23] 
we will discuss in more detail the possible mechanisms of he- 
lium and hydrogen burning layers and expulsion efficiencies. As 
the mass of the WD approaches the critical mass limit of MsNia = 
carbon ignites. Due to the degenerate nature of the WD 
a thermonuclear runaway takes place, leading to a SNIa event. 

2.2. Common envelope phase 

The classic al way to parametrise the CE p hase is by the a- 
formalism (iPaczvnskill 19761 IWebbinklll984t) which is based on 
conservation of energy: 
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where the left hand side represents the binding energy of the CE 
Aiibind and the right hand side a fraction a of the orbital energy 
AEorb of the binary. M represents the mass of the donor star and 
the subscripts e and c refer to the envelope and core of the donor 
star. The secondary mass is denoted by m and is assumed not to 
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Fig. 1. Binary evolution for a single-degenerate SNIa progenitor 
In this case the donor is a MS star. The top and bottom parts 
of the figure have different scales due to the common envelope 
phase, denoted as CE in the figure. 

vary during the CE phase. The a, and a/ are the initial and final 
separation of the binary. R is the radius of the donor star and A 
is a structural parameter of the envelope. All the uncertainty can 
be taken into one parameter, the product a A. For higher values of 
aA, the CE is expelled more efficiently and the change in orbital 
separation will be less dramatic. However, in all cases the orbital 
separation shrinks as a result of the CE phase. Contradicting this, 
obse rved double WD binarie s seem to have mass ratios close to 
one dMaxted & Marsh 1999). To form a double WD pair the bi- 
nary usually goes through two CE phases. One when the primary 
star evolves into a giant and subsequently a WD and one when 
the secondary evolves. In order for the mass of the WDs to be 
roughly equal, the conditions at the start of the first CE cannot 
differ much from those at the onset of the second CE. Therefore 
the orbital separation of these systems cannot shrink substan- 
tially during the first CE phase. 

As a possible solut ion to this problem the y-algorithm was 
introduced by Nelema ns et al.l (1200 0). Based on the conservation 
of angular momentum, angular momentum is lost in a linear way 
with mass according to: 
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where J represents the angular momentum of the binary, M,o, 
is the total mass of the binary and the other parameters have 
the same meaning as before. Depending on the mass ratio, the 
orbital separation will be unchanged or decreased due to the de- 
scription of the CE phase. The physical mechanism behind the 



y-des cription remains unc lear Interesting to note here is that re- 
cently I^^Ddretal] (1201 2h suggested a new evolutionary model 
to create DWDs. These authors find that mass transfer between a 
red giant and a MS star can be stable and non-conservative. The 
effect on the orbit is a modest widening, with a result alike to the 
y-description. 

2.3. Growth of the white dwarf 

How efficiently the WD grows in mass as a result of the accretion 
is described by the total retention efficiency ?/tot- This represents 
the fraction of transferred hydrogen-rich matter from the c om- 
panion that eventually burns into carbon-rich matter and stays on 
the WD. Therefore Mwd - '7tot|A^fcompl, where Mwd represents 
the mass growth rate of the WD and Mcomp the mass transfer 
rate of the companion. Note that Mcomp < 0, while Mwd > 0. 
Ultimately the retention efficiency of the transferred mass deter- 
mines in which systems a SNIa event takes place and therefore 
shapes the SNIa rate and delay time distribution. 

2.3.1. Hydrogen-rich accretion 

When hydrogen-rich material is transferred onto a WD, the evo- 
lution of the accreted layer depends on the mass transfer rate 
of the companion Mcomp and the WD mass Mwd- The evo- 
lution can be split into three distinct regimes, separated by 
a so-called steady and crit i cal accretion rate Mst (MwD) and 
Mcr(MwD) (lNomotdll982bt iHachisu & Katdl200lh . These are 
both functions of the WD mass and are discussed in detail be- 
low. 

For low mass transfer rates (IMcompl < M^t) the hydrogen- 
rich matter is transferred to the WD conservatively, so that the 
mass transfer rate onto the WD is Mh = IMcompl- The tempera- 
ture in the accumulating hydrogen layer on the WD is too low to 
ignite the hydrogen immediately. The layer will gradually grow 
in mass, increasing the pressure and temperature until hydrogen 
ignition values are reached. Then nuclear burning will quickly 
spread through the layer , which starts expand ing. The strength 
of these so-called novae ( IStarrfield et al.|[l972l) depends strongly 
on the amount of mass involved. This also determines how much 
of the processed layer is expelled from the binary and how much 
falls back onto the WD as helium. The accretion rate of the 
helium-rich matter will be denoted by Mhc- We define the hy- 
drogen retention efficiency to be the fraction of transferred mass 
which is burned and retained by the WD: 



Ma 
Mh 



(3) 



If none of the matter stays on the WD rjn - and if all of the 
transferred matter is retained by the WD as helium t^h = 1 ■ In the 
nova regime ?7h < 1 • For very strong novae not only the complete 
layer but also some of the WD itself may be expelled. In these 
cases T]H < 0. 

Intermediate mass transfer rates are those between the steady 
and critical value, so that Mjt < IMcompI ^ M^ . All of the mass 
lost by the companion star is transferred to the WD, so that Mh = 
l^compI- On the WD the hydrogen accumulates in such a way that 
the ignition values are reached at the bottom of the layer and the 
hydrogen stably bums into helium. No material is lost from the 
layer or the WD, therefore ?7h = 1 and Mhc - Mh. 

For higher mass transfer rates exceeding the critical value 



> M„) a hydrogen red-giant-like envelope forms 



around the WD. If the high accretion rate is maintained a CE 
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phase will soon follow. However, hydrogen burning on top 
of th e WD is strong enough for a wind to develop from the 
WD dKato & Hachisu 1994^HKN99b). This attenuates the mass 
transfer rate Mh and hydrogen retention efficiency rju, and en- 
sures that the CE phase is postponed. Part of the hydrogen-rich 
matter, corresponding to M^-, is burned into helium and stays on 
the WD. The remainder Mh - Ma- = M„ind is carried away by the 
wind. In this wind regime, Mue = M„ = rjuMu and so /ye < 1 . 

A possible second effect influencing the retention efficiency 
in this high mass transfer rate regime, results from the wind if 
it reaches the companion star There it heats up the envelope of 
the companion, which expands and can then be stripped off by 
subsequent wind (HKN99b). It is now no longer the case that all 
of the mass lost by the companion is transferred to the WD. The 
part that is transferred to the WD is equivalent to the amount lost 
by the companion minus the amount stripped away by the WD 
wind: Mh - IMcompI - lA^^smpI, where Msmp < 0. The strength of 
this stripping effect is represented by the stripping parameter ci , 
defined by 



Mstrip = CiMwind- 



(4) 



Details on how and when this stripping effect is taken into 
account are described in section|3]and appendixlAl 

2.3.2. Helium-rich accretion 

After a fraction of the hydrogen-rich matter is burned into 
helium-rich material, the helium-rich material accumulates in a 
layer on the surface of the WD. This layer sits in between the 
WD and the (still accumulating) hydrogen layer For the helium- 
rich accretion similar regimes exist for low and medium he- 
lium accretion rates, resulting in helium novae and steady he- 
lium burning respectively. No helium wind regime exists be- 
cause these high helium accretion rates cannot be reached by 
the hydrogen burning. The fraction of helium that is burned into 
carbon-rich material is represented by the helium retention effi- 
ciency rjue- The total retention efficiency is then given by; 



3. Method 



(5) 



To test the viability of SD SNIa theory a population of binary 
stars is simulated and the corresponding theoretical SNIa rate is 
compared to observational re sults. We employ the binary pop- 
ulation synthesis code SeBa ( Portegies Zwart & VerbuntI [19961 
iNelemans et alJl200lblToonen et alJl2012h for fast stellar and bi- 
nary evolution computations. In SeBa stars are evolved from the 
ZAMS until remnant formation. At every timestep, processes 
as stellar winds, mass transfer, angular momentum loss, mag- 
netic braking and gravitational radiation are taken into account 
with appropriate recipes. SeBa is a Monte Carlo based code in 
which the initial binary parameters are generated randomly ac- 
cording to appropriate distribution functions. The initial mass of 
the primary stars is drawn from a Kroupa Initial Mass Function 
(iKroupa etal.lll993h which ranges from O.I-IOOM0 and the ini- 
tial mass of the secondary from a flat mass ratio distribution with 
values between and 1 . We assume that a complete stellar pop- 
ulation consists of 50% binary stars and 50% single stars. The 
semi-major axis of the binary is drawn from a power law dis- 
tribution with an exponent of -1 (lAbtll 19831) . ranging from to 
10^ Rq and t he eccentricity from a thermal distribution, ranging 
from to I (lHeggidlI975h . For the metallicity solar values are 
assumed. 



In this research we differentiate between four ways to deter- 
mine which binaries will give rise to a SNIa event. For the first 
method (model NSKH07) we adopt the retention effic i encies 
given by NSKH07, including updates from lHachisu et all ( l2008l 
hereafter HKN08). Using SeBa in combination with these reten- 
tion efficiencies, a population of binaries is simulated and those 
systems in which a C/O WD reaches the critical mass limit are 
selected as SNIa progenitors. For the second appr oach (mo del 
RBF09) we utilise the retention efficiencies used in IRBF09I and 
for the third approach (model YIO) the efficiencies based on 
YIO, The assumed retention efficiencies differ strongly. For the 
final method (model Islands) the phase of stable mass transfer 
onto the WD is not modelled by SeBa. Instead we select those 
systems that at the formation of the C/O WD lie in specific re- 
gions of the parameter space of orbital period, WD mass and 
companion mass (Porb - A^wd - A^comp)- Systems in these spe- 
cific regions, hereafter islands, are determined to lead to SNela 
(HKN6s : ,HKN99ta . and references therein). 

In order to compare our simulated rates to the results of re- 
search groups that use the same retention efficiencies, the pre- 
scription of the CE is varied to match different BPS codes. We 
distinguish between two models. The first model is based on 
the cy-prescription assuming a value of a - I and A = 0.5 as 
in accordance with lRBF09l The secon d model assumes the y- 
formalism with a value of y = 1.75 (see INelemans et a"Dl200lh FI 

3.1. Retention efficiencies 




M(Mo/yr) 



Fig. 2. Examples of hydrogen and helium retention efficiencies 
as a function of the mass transfer rate for a WD of IMq. The 
prescriptions for t he hydrog en retention efficiency r/ u] (sohd 
Une - INSKH07I : iHKNOSi) & i]m (dashed line - IrBF09I : 
YIO|) and the helium retention efficiency ?7H ei (dot-dashed line - 
Hachisu et al."l999a':' Kato & Hachisulll999h & J]He2 (dotted line 
- Iben & Tutukov 19961) are shown in AppendixlAl Note that the 
position of the peak of rjui is dependent on the exact prescription 
for Mci- The mass transfer rate for hydrogen accretion is given 
by IMcompI- For helium accretion this is Mee- 



In this section we describe the three prescriptions for the hy- 
drogen and helium retention efficiencies as given by INSKH07I 



' The y-formalism is applied, unless the binary c ontains a com pact 
object or the CE is triggered by a tidal instability (iDarwi nl lTST^ for 
which the a-formalism is used. 
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IRBF09I and lYlOl The total retention efficiency is a combina- 
tion of the hydrogen and heHum one as described by eq.|5] The 
retention efliciency used in NSKH07, including updates from 
iHKNOa is a combination of the hydrogen retention efficiency 
rjm and the helium retention efficiency zyHei, see Fig.|2] For 77H1 
the region of steady burning occurs at mass transfer rates be- 
tween 



Ms, = 3.1 ■ 10"^(MwD - 0.54) Moyr"' and 



7.5 



10-'(MwD 



• 0.40) Moyr" 



(6) 



(7) 



which are the updated formulae from 'HKN08I with MwD in units 
of Mq. This region is visible in Fig.|2] as the plateau at inter- 
mediate mass transfer rates. In the wind regime |Mcompl > M„, 
the wind-stripping effect is included wit h ci = 3 (see eq.Bl). For 
the prescr iption of ?7Hei see Appendix lAldKato & Hachisul 19991 : 
lHKN99bh . 

The retention efficiency based on lRBF09l combines t/hc 1 with 
77h2, see Fig.|2] In this model there is no region of steady hy- 
drogen burning. Instead the nova regime immediately borders 
the wind regime. The critical mass transfer rate is given 
by eq.|7] The stripping effect is not taken into account, hence 
ci = (see eq.|4]i. The retention efficiency in the nova regime 



compl 



< Mcr) is base d on an interpolation of results from 



iPrialnik & Kovetzl (ll995h . Our own fit to these results is included 
in AppendixlAl For lower (higher) Mwd than assumed in Fig.|2] 
the retention efficiency in the nova regime stays the same be- 
tween 10"** < IMcompI < 10"^ and shifts with the peak in the 
regime 10"^ < |Mcomp| < M„. 

The retention efficiency based on lYlOl combines rim with 
J]He2, see Fig.|2] The critical rate for the hydrogen accretion is 
given by 



M„ = 10 



-9.31+4.12MwD-l-42M^, 



™ Moyr 



with Mwd in units of Mq. For the retention efficiency 
in the hydrogen nova regime the same interpolati on from 
IPrialnik & Kovetzl (Il995h is used as for the efliciency of lRBF09L 
see AppendixlAl Again there is no region of steady hydrogen 
burning, the wind is taken into account but the stripping effect is 
not. Due to the different prescription for M„ the position of the 
peak is slightly different from the example curve. Using eq.|8]as 
the border to the wind regime for hydrogen accretion, the peak 
occurs at smaller values of IMcompI for a given Mwd, as compared 
to eq.|2l This reduces the total retention efficie ncy ?7tot substan- 
tially, as we will see in Fig. [3] Following YIO, the prescription 
for 77He2 was taken from Jben & T utukov (1996). 

The exact presc riptions of the hyd roge n and helium reten- 
tion efficiencies for INSKHOI IrBFO^ and IyToI are detailed in 
AppendixlAl Fig. [3] shows the large variety in the total retention 
efficiency for the three prescriptions as a function of mass trans- 
fer rate of the companion . The most optimistic retention effi- 
ciency is that of lNSKHOTl and the most pessimistic that of lYlOl . 
Note that in the regime with strong novae the retention efficiency 
T/tot ^ 0. Here we have set T/tot = in this regime for simplicity 
as these systems are not part of the SD channel. 

A complication of this method is that the instantaneous mass 
transfer rates in binary population synthesis codes, such as SeBa, 
is only an approximation to the one derived from detailed stel- 
lar evolution codes. Therefore, we also implemented a hybrid 
method in which progenitors are selected from binary popula- 
tion synthesis according to results from the literature based on 
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Fig. 3. Total retention efficiencies, resulting from different com- 
binations of the hydrogen and helium retention efficiencies. In 
this figure we have assumed Mwd = IM© as an example. The 
mass transfer rate M = IMcompI- For lower (higher) Mwd the 
maximum retention efficiency shift to lower (higher) IMcompI- 
For a more detailed explanation of the retention efficiencies see 
Sect.lO 



more detailed mass transfer tracks. These tracks can be cal- 
culated by an analytical approac h e.g. iHKN99bi or by a de - 
tailed binary evolution code e.g. iLi & van den Heuvell (l2002h . 
iHan & Podsiadlowskil (|2004 and lHKNOSl . 



3.2. Islands 



(8) In this method, the progenitor binaries are selected on P, 



Mwd, and Mcomp at the time of WD birth (third stage in Fig.[T}. 
The parameter regions used are shown at the right hand side in 
Fig-Sfor binaries containing a WD and a MS star or Hertzsprung 
gap star (WD+MS binaries) and Fig.jS] for those containing a 
WD and a red giant star (WD +R G binaries ). They are made to 
match the islands in IHKN08I and lHKN99bl respectively, which 
are shown on the left hand side in these figures. For the WD-hMS 
channels a moderate wind-stripping parameter of ci =3 is taken 
(see eq.|4]i. 



WD + MS system (cl = .3) 

M^[s.-0.7. 0.8. 0.9. 1.0. 1 1 M„ 




WD+MS: Mwd = 0.7, O.B, 0.9, 10, 1.1 M. 




Fig. 4. Initial parameter regions for the WDh-MS track to SNela. 
The different contours are for different WD masses (increasing 
in size wit h increasi ng mass) and the stripping parameter ci = 3. 
Left: from iHKNOa The hatched regions indicate SNIa explo- 
sions with short delay times of f < lOOMyr for M = Q.IMq and 
M = 1 . 1 Mq. Right: initial parameter regions as used in this work. 
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WD+RG:M„oo^0.a. 0.9, 1.0,1.1 ^^| 




WD+RG: MwD = 0.8, 0.9, 1.0, 1.1 M,: 




Fig. 5. Initial parameter regions for the WD +RG track to SNela. 
The different contours are for different WD masses (increasing 
in size wit h increasin g mass) and the stripping parameter ci - I. 
Left: from lHKN99bl Right: initial parameter regions as used in 
this work. 



4. Results: SNIa rates and delay times 

We have simulated a population of SD SNIa progenitors for the 
four approaches described in Sect.|3]assuming our model for the 
y-prescription or the a-prescription. The delay time distributions 
for the y-prescription simulations can be seen in Fig. |6] and for 
the a-prescription in Fig.|7] The delay time t is the time at which 
a SNIa occurs, where f = is the time when the binary is born 
as a double ZAMS-star binary. A DTD shows the distribution of 
delay times after a single starburst involving a large number of 
binaries. For model Islands, we assume the phase between the 
onset of mass transfer and the SNIa explosion is short compared 

to the lifetime of the binary. 

The progenitor islands are based on the work of INSKH07I 
andiHKNOS, as is our model NSKH07 that uses the retention ef- 
ficiencies directly. If we compare the DTD resulting from model 
Islands with model NSKH07, there is a noticeable difference at 
early delay times. The DTD of model Islands peaks earlier and 
higher. The reason is the extend of the Islands to high donor 
masses of about 6M0for WD+MS progenitors. These massive 
companions fill their Roche lobes soon and therefore the SNIa 
explosion occurs at earlier delay times (f < IGyr). In SeBa, bi- 
naries with such high mass ratios undergo unstable mass transfer 
and a CE phase. They do not develop a SNIa explosion, resulting 
in less SNela at early delay times from the retention efficiencies 
method. 

Our progenitor regions in model Island consist of adja- 
cent slices that increase in size for increasing Mwd, for both 
WD+MS and WD+RG progenitors. These slices have a thick- 
ness of O.IMq around the values Mwd = 0.7, 0.8, 0.9, 1.0 
and l.lMo for WD+MS binaries and Mwd = 0.8, 0.9, 1.0 and 
I.IM0 for WD+RG binaries. Ideally the volume enclosing the 
progenitor binaries has smooth edges in all three dimensions. 
Since this is not the case in our approach we explored the re- 
sulting inaccuracy by increasing our resolution. We double the 
amount of islands covering the same range of WD mass, so the 
thickness of each island slice is now O.OSM©. We found that the 
change in the final integrated SNIa rate is 3-4% and the DTD is 
barely affected. 

The SNIa DTD and integrated rate are strongly influenced 
by the prescription of the retention efficiency, as shown in Fig.|6] 
and Fig.|2l and Table[T] For more pessimistic values of the to- 
tal retention efficiency the rates and overall height of the DTD 
decreases. Note that the total retention efficiency of model YIO 
is so small that no SNela developed and only an upper limit is 
given in Table[T] Model NSKH07 and RBF09 give rise to DTDs 
that differ most strongly at short delay times. This is due to high 
donor masses that transfer mass at a high rates, where the total 




8.0 8.5 9.0 

Delay time (yr) 



10.0 



Fig. 6. Delay time distribution assuming the y-algorithm for 
the CE with y - 1.75. Different lines correspond to the delay 
time distributions that result from four different approaches de- 
scribed in Sect|3] all for the y-prescription. The dot-dashed line 
shows the result of the Islands selection method, the three other 
linestyles correspond to the retention ef ficienc ies as in Fig. [3] 
Note that for the retention efficiency of lYlOl no SD binaries 
evolve into a SNela. 



retention efficiencies differ significantly, see Fig .[3l This re sults 
from the inclusion of the wind stripping effect in lNSKHOTl . 

Changing the prescription used for the CE phase modifies 
the integrated SNIa rate by a factor of about two. When the 
y-prescription is applied, less (very) close binaries are created 
as the binary orbits do not shrink as effectively as for the a- 
prescription. An exception to this is when the mass ratio of the 
binary is extreme at the moment of Roche lobe overflow. These 
systems merge at short delay times t<0. IGyr. Therefore the DTD 
when the y-prescription is assumed declines faster with time 
than when the or-prescription is assumed. The classical evolu- 
tion path towards a SD binary as depicted in Fig.lT] is less com- 
mon in the BPS model using the y-prescription. Assuming the a- 
prescription 80% of the SD binaries evolve through this channel, 
where as for the y-prescription this has decreased to 30-40% (de- 
pending on the retention efficiency). In the most common evo- 
lution path for the model assuming the y-prescription, the first 
phase of mass transfer is stable. 



5. Discussion 

The theoretical SD SNIa rates that follow from SeBa in this 
study can be compare d to the resul t s of various othe r BPS 

iMennekens et all I2OIOI: 



research groups (e.g. IRBF09 



IWang et al.ll20IOt IClaevs et al.ll20I Ih . Tabled shows the time- 
integra ted SNIa rates of these groups taken from lNelemans et all 
(l20I2h . The disagreement in the rates is large up to a factor 
of about 600, but so far no explanation has been found. Note 
that Nelemans et al. (20121) rescaled the results (if needed) to the 
same initial distribution of parameters as discussed in Sect. [3] 
Assumptions and simplifications vary between the different BPS 
codes (See Toonen et al. in prep, for a study on this) causing 
differences in their predictions. An important assumption for 
the BPS simulations is the assumed CE-prescription, which is 
also given in in Table|2] The effect of different values for a 
has been studied by several groups, e.g. the effect on the DTD 
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Table 1. Time-integrated SNIa rates in the SD channel for different combinations of the retention efficiency and the common 
envelope prescription in units of 10"^Mo~'. 



Model 


Approach used 


y-prescription 


ff-prescription 


Model NSKH07 


Retention efficiency of NSKH07' 


0.59 


1.3 


Model RBF09 


Retention efficiency of RBF09- 


0.19 


0.35 


Model YIO 


Retention efficiency of YIO^^ 


<0.001 


<0.001 


Model Islands 


Islands of HKNOS"*, HKN99b'' 


0.73 


1.5 



Notes. ' iNomotoetall ( 120071) : ^ iRuiter et al] ( 120091) : ^ lYungelsonl j201(J ): ^ iHachisu et all f200l) : ^ iHachisu et all <1999br) 




8.0 8.5 9.0 

Delay time (yr) 



10.0 



Fig. 7. Same as in Fig. |6] but assuming the o'-formalism with 
a A - 0.5. Again no SD b inaries evolve into a SNela for the 
retention efficiency of lYlOl 



dWang et al. 



IClaevs et al 



20101: iRuiter et alj [20091: iMennekens etal] 120101: 



201 ih . The effect on the integra ted rate of a small 



change in a is of the o rder of a factor 0.7-3 (iRuiter et al.ll2009l 
IMennekens et alj2010l) . however the effect can be up to an order 
of magnitude for larger changes in a (Claeys et al. in prep.). 

The entries in Table|2] are ordered to increase in rate. The 
smallest rate is from the work of Yungelson (IyToI) . Similarly we 
find the lowest rate for model YIO, however we find an even 
lower rate than Yungelson. The preliminary rate of IClaevs et alj 
(1201 lb i s significantly lowe r than our corresponding model. 
However lClaevs et al.l (1201 ih use a variable lambda prescription, 
whereas in our bes t correspond i ng mo del we set aA - 0.5. The 
integrated rates of IRuiter et al.l (l2009b is a factor of about two 
lo wer than the rates of our best corresponding mod el. The rates 
of IMennekens et al.1 (1201 Ol) and IWang et~al] j2010l) are a factor 
of a about seven and two higher respectively. So although our 
results do recover the trend in the theoretical predictions from 
different binary population synthesis codes, they do not fully ex- 
plain the large disagreement among them. 

The integrated rates based on observations are given in 
Table|2] The most recent measurements are based on field galax- 
ies and generally show lower rates, while earlier estimates based 
on galaxy clusters are higher At this moment it is unclear if 
the different observed integrated rates are due to systematic ef- 
fects or if the re is a r eal enhancement of SNela in cluster galax- 
ies. See also Maoz et al . (2012) for a discussion. Even though 
the different retention efficiency models affect the SNIa rates 
with a factor > 10^, none of the integrated rates comfortably 



matches the observed rates, especially those from galaxy clus- 
ters. In addition, the DTD reconstructed from observations typ- 
ically show a continuation to longer delay times, which are ab- 
sent in all our SD DTDs. We conclude from this, that in the 
current model of SD SNIa theory, the main contribution to the 
SNIa rate comes from other evolution channels. One possible 
channel involves semi-detached binaries in which a WD ac- 
cretes fiTOi]_a_Jwdrogei>poo£jiel^^ such as sdB 



stars. 



Wang et aU ([2009al) . IRuiter et al.l (|2009|) and IClaevs et~al1 



(201 II) showed that in this channel the DTD peaks at delay times 
of about lOOMyr, although rates at this delay time vary between 
10""* - lO"-yr"'(lO'°M0)"'. The contribution from binaries of 
the DD channel is debated heavily. Explosion models favour 
accretion-induced colla pse to a neutron star over a SNIa event 
(see e.g. .Hillebrandt & Niemevej|2000t for a review). However, 
BPS codes find more SN Ia events from the standard DD chan- 
nel than the S D chann el (IrBFTOI: I YlOl: IMennekens et alll2010l: 
iToonen et al ]|2012fl T oonen et al. I (l2012h study the contribution 
from the double-degenerate channel with SeBa comparing the 
a- and y-prescription for the CE. They find that even though 
the DD DTDs fit the observed DTD beautifully, the normalisa- 
tion does not by a factor of about 7-12 compared to the clus- 
ter rates. Taking into account the new rates from field galax- 
ies, the factor becomes about 1.2-12. Other contributions to the 
SNIa rate can possibly come from e.g. core-degenerate mergers 
(iKashi & Sokeil 1201 Ij), double detonat ing sub-Chandrasekhar 
accretors (see e.g. Kromer et al. '2010") or Kozai oscillations 
in tri ple systems tShappee & Thompsoft ,2012; .Hamers et aTl 
120131) . 



6. Conclusions 

In this work we have studied the effect of the poorly under- 
stood phase of WD accretion in the context of supernova type 
la rates. We employed the binary population synthesis code 



SeBa dPortegies Zwart & VerbuntI [19961 iNelemans et all 120011: 
iToonen et al.l 2012 ) to study the SNIa rates and progenitors 



for different CE prescriptions. We differentiated between four 
models a s suming eith er on e of three retention efficiencies of 
INSKH07L[RBFTOI andFYTol or making a selection of SNIa pro- 
genitors based on binary parameters at the time of WD formation 
(HKN08; HKN99b). The three retention efficiencies assumed by 
different binary population synthesis codes differ strongly. The 
difference comes from a lack of understanding at low mass trans- 
fer rates where novae occur, and mass transfer rates higher than 
the rate for stable burning. This is true for the accretion of hydro- 
gen that is transferred to the WD by the companion, as well as 



^ Note that the int egrated rat e from the violen t merger model for dou- 
ble WDs ( Pakmor e tall lIoTol [2 011. 2012; Ro pke et alJiMl ) can be 



- ,2h. 
much lower tToonen'et al.l2012l : IChen et al.ii2012LlRuiter et al.l2012h 
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Table 2. Time-integrated SNIa rates in the SD channel for the assumptions of different research groups in units of 10 "^Mq"' 



BPS 


WD accretion 


SNIa rate 


CE 


SNIa rate 


CE 


Yungelson' 


YIO 


<0.001 


7 = 1.5 


0.006 


7= 1.5 


Claeys^ 


HKN99b 


1.3 


aA = 0.5 


0.13 


a = I, A = variable 


Ruiter* 


RBF09 


0.35 


aA = 0.5 


0.17 


aA = 0.5 


Wang/Han" 


HKN99b 


1.3 


aA = 0.5 


2.8 


aA = 0.5 


Mennekens'' 


HKN99b, HKN08 


0.55 


aA = 1 


3.7 


aA = 1 


Observed^ 




4-26 







Notes. Columns 3 and 5 shows the integrated rate predicted by SeBa and the BPS group in question respectively. The a ssum pti on of the BPS group 
for the CE-evolution is shown in column 6 and the be st corresponding CE-mo de l in SeBa is shown in column 4. ' | Y10|: ^ Claevs et al. ! (1201 ih 
(preliminai-y results); ^ IRBFO^ : Imng et all ( l20loh : ^ iMennekens et al.l ( l201(ih : ^ iMaoz & Mannuccil ^201 Ih ^ lPerrett & et al.l7201Z) : ,Maoz et all 
(l2012h : lGraur & Maozld2'oi2h 



the accretion of helium that has been burned from hydrogen on 
the WD. The efficiency with which a C/O WD grows in mass is 
strongly affected by the combination of the efficiency of hydro- 
gen and helium accretio n. For example the hydrogen and helium 
retention efficiencies o f lYlol are maximal at different ranges of 
the mass transfer rate resulting in a low total retention efficiency. 

The total number of SNIa progenitors is significantly influ- 
enced by the choice of the model. The integrated SNIa rate vary 
between < 1 ■ 10"^ - 1.5 • 10"'' M©"', where the rates a re highest 
for the model that assume the retention effici encies of lNSKHOTI 
and only an upper limit can be given forlYlol Our method based 
on the pa rameter space of binaries at WD birth of |HKN0"8| and 
lHKN99bl consists of a discrete set of islands in (Poib - M^f^o - 
Mcomp)- The discretisation introduces an error of 3-4% on the in- 
tegrated rates. We showed in this paper that independent of the 
WD accretion model the SNIa rate approximately doubles when 
the ff-prescription is assumed with aA = 0.5 as compared to the 
y-formalism with y = 1.75. The effect of different values for 
a on the SD SNIa rate is a factor of 0.7-3 (iRuiteret al.ll2009t 
IMennekens et alJbOlOh for small changes in the CE-efficiency 
and up to an order of magnitude for larger changes in the CE- 
efficiency (Claeys et al. in prep.). Also note that throughout this 
work we have assumed solar values for the stellar metallicities. 
Using a broad range of metallicities in the BPS code might influ- 
ence the SNIa rates and DTDs, both in th e Island s and retention 
efficiency approach. See lKobavashi et al.l(ll998bl) for a study on 
how the progenitor islands depend on metallicity. 

Several predictions for the SD rate have been made by dif- 
ferent groups using binary population synthesis simulations. The 
results show a wide spread and do not agree with each other 
or with observational data. The integrated rates vary between 
6 • 10"^ - 3.7 • lQ-'^M0-\ In this study we find that also the 
model for the WD accretion is a major source of uncertainty on 
the SNIa rates. The different prescriptions for the retention ef- 
ficiency introduce a source of uncertainty with an effect on the 
integrated rates by a factor of about 3-4 (comparin g ou r models 
that assume the retention efficiencies of lNSKHOTi with lRBF09h 
or even larger u p to a factor of a few hundred (comparing with 
the efficiency of lYlOh . Although our results do recover the trend 
in the theoretical predictions from different binary population 
synthesis codes, they do not fully explain the large disagreement 
among them. As the exact origin of the differences in the SD 
rate remains unclear, Toonen et al. in prep, study the results of 
four different binary population synthesis codes and investigate 
the importance of the different assumptions and numerical ap- 
proaches in these codes. 
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Appendix A: Retention efficiencies 

The total retention efficiency is the product of the hydrogen and 
helium retention efficiencies. In the equations in this appendix, 
aU M are in units of Mo/yr 

A.1. Retention efficiencies based on \NSKHO^ 

The hydrogen retention efficiency is a strong function of the 
mass transfer rate of the companion Mcomp- Three regimes can be 
distinguished, separated by the stable and critical mass transfer 
rates [M^t and Mcr), see Table lA. Il In the nova and stable burning 
regime all the mass lost by the companion is transferred onto the 
WD. In the stable regime all hydrogen-rich matter is burned into 
heUum-rich matter and all stays on the WD, so that rjn = I. The 
nova regime is linearly interpolated between the lower bound- 
ary at IMcompI = 10"^ and the start of the stable regime at M^t- 
In the third regime the nuclear hydrogen burning on top of the 
white dwarf is so strong that a wind is produced which not only 
attenuates the mass transfer rate but can also strip the compan- 
ion's outer envelope. It is no longer the case that all the mass lost 
by the companion accretes onto the white dwarf. The maximum 
that the white dwarf can accrete is M„, which is a fraction rju 
of the mass transferred to the white dwarf Mh- All the excess 



is blown off by the wind. The amount of matter that is stripped 
from the companion is defined by eq.|4]and the stripping param- 
eter is taken to be ci =3. 



Table A.l. The three reg imes for di fl'erent mass transfer rates of 
the companion star, as in lNSKHOTi 



Mcomp range 


Nova regime 


IMcompI < M^t 


Stable burning regime 




Wind and stripping regime 


IMcompI > M^T 



The exact prescriptions for the stable and critical mass trans- 
fer rates are: 



Mst = 3.1 • 10"^(^^^ - 0-54) and 



Mc- = 7.5 ■ 10 



(A.l) 



(A.2) 



For the hydrogen retention efficiency we arrive at the follow- 



ing: 



(log(IMcompl) 



if IMcomp 



< 10-^ 
7)/((log(M,0 + 7) 

if 10-^ < IMcompl < M„ 
if Mst < IMcompl < M„ 

compl 



(A.3) 



McVMh if Mer < IMeomol < 10-4 



with, for the wind and stripping regime, 

(ci + DM,, 



Mh ciM„- + |M, 



(A.4) 



compl 



This last equation follows from some algebra, taking into 
account that Mcomp, M„ind and Mstrip < because they describe 
matter travelling away from one of the stars and using eq.S 

For the helium retention efficiency the following prescrip- 
tions were used: 



if Mhc < lO-^ '* 
-0.175(log(MHe) + 5.35)2 + ^ Q5 

T/He = i if lO-^ ** < Mhc < lO-^ "* 

1 if 10-^ '^ < Mee < lO"-*' " 
if Mhc > 10-5 ° 



(A.5) 



where Mhc - tJuM, 



comp- 



A.2. Retention efficiencies based on \RBFO§i 

For the hydrogen retention efficiency two regimes can be distin- 
guished, the nova regime and the wind regime, see Table lA. 2l 
The stripping effect is not taken into account (ci = 0). In 
the hydrogen nova regim e an interpolation of the results from 
iPrialnik & Kovetzl (Il995h is used for t/h. The amount of mass 
transferred to the white dwarf is always equal to the amount of 
mass lost by the companion, |Mh| - IMcompl- 
The critical mass transfer rate is given by: 



Mc, 



: 7.5- 10 



1 



M» 



•0.40 



(A.6) 
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Table A.2. The two regi mes for different mass transfer rates of with 
the companion star, as in lRBF09l 



Nova regime IM^orapl < M„ 
Wind regime IMcompI > ^cr 



Me, 



IMcompI 

The helium retention efficiency is: 

if Mee < 10"^-5 



The hydrogen retention efficiency is: 



if IMeompI < 10-8 

0.25(log(|Meomp|) + 8) if 10-« < |Meo,npl < 10 



?7He = 



^ if 10-^-5 < Mhc < 10- 
0.95 if 10-5-^ < Mhc < lo- 



ur- 



if 10 



-5.4 



He 

< Mee < 10- 



0.25 + 0.75(log(|M,omp|) + 7)/(log(M„) + 7) 

if 10-^ < |Meo,npl < 

M„/Mh if M„ < IMeompI < 10-4 



if MHe > 10 

here Mhc = 'ynMcomp- 



-4.0 



where 

M„ M„ 



Mh |M,on,p| 

The helium retention efficiency is given by 

if Mhc < 10-^-3 

-0.175(log(MHe) + 5.35)2 + j 
/7He = \ if 10-'-^ < MHe < 10-^" 

1 if 10-5 '' < Mne < 10-5 ° 
if Mee > 10-5 ° 



(A.8) 



(A.9) 



where Mhc = ?7HMc omD. Note the similarity with the helium re- 
tention efficiency of lNSKHOTl except for the lower limit of the 
mass transfer rate of the stable burning regime. 



A.3. Retention efficiencies based on \Y1(X 

Two regimes can be distinguished for the hydrogen retention ef- 
ficiency, the nova and wind regimes, see T able lA.3l The same 
interpolation from lPrialnik & Kovetzl (1 19951) is used for the hy- 
drogen nova regime. The stripping effect is not taken into ac- 
count (ci = 0). In all cases the amount of mass transferred to the 
white dwarf is equal to the amount of mass lost by the compan- 
ion, |MhI - IMcompI- It is very similar to the retention efficiencies 
in Sect. lA.2l except that the prescription for M„ is different. 



Table A.3. The two regi mes f or different mass transfer rates of 
the companion star, as in lYlOl 



Mcomp range 

Nova regime IM^orapl < M^, 
Wind regime IMcompI > Mc, 



The prescription for the critical mass transfer rate is: 

M„ = 10-''-31+412MwD-1.42M^„^ jQ) 

The hydrogen retention efficiency is given by: 

if IMeompI < 10-8 

0.25(log(|Meompl) + 8) if 10-8 < IMcompI < 10-^ 

7?H = i 0.25 + 0.75(log(|Mcompl) + 7)/(l0g(Mcr) + 7) (A. 11) 

\f 1 r\-7 ^ \A/f I 

compi 
■ comp I 



if 10-^ < IMcompI < Mer 

McJMh if Mcl- < IMcompI < 10-4 
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